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Fractional transport of particles on a comb structure in the presence of an inhomogeneous con- 
vection flow is studied. The large scale asymptotics is considered. It is shown that a contaminant 
spreads superdiffusively in the direction opposite to the convection flow. Conditions for the realiza- 
tion of this new effect is discussed in detail. 



PACS numbers: 05.40.-a, 05.40.Fb 

A realization of superdiffusion on a subdiffusive me- 
dia, e.g. on a comb structure is an example of the 
fractional transport due to the inhomogeneous convec- 
tion. These studies of the space-time evolution of an 
initial profile of particles in a specific media due to the 
inhomogeneous convection flow are arisen in a variety of 
applications such as transport of external species (pollu- 
tion) in water flows through porous geological formations 
H> m ! Problems of diffusion and reactions in porous cat- 
alysts Q and fractal physiology The conditions for 
the inhomogeneous convection which is responsible for 
the superdiffusive spreading of an initial packet of parti- 
cles have been found in Q. A classification of possible 
scenarios of the space-time evolution of a contaminant 
in the presence of the inhomogeneous convection on the 
comb structure is presented there. The external forcing 
has been taken in the power-law form for the convec- 
tion current j x {t,x,y) — vx s S(y)G(t,x 1 y), where a dis- 
tribution G(t, x, y) describes the evolution of the initial 
profile, while vx s S(y) is the inhomogeneous convection 
velocity. When s < subdiffusion is observed 0, Q- 
When s > it is superdiffusion. The homogeneous con- 
vection with s = corresponds formally to the normal 
diffusion, but the effective diffusion coefficient is deter- 
mined by the external forcing v 0. The frontier case 
with s = 1 corresponds to the log-normal distribution of 
transport particles, where one deals with not a sum of in- 
dependent random variables but with their multiplication 
When s > 1 the fractional transport corresponding 
to superdiffusion possesses specific features. This case is 
the main subjective for the investigation in the present 
Letter. We consider both the fractional transport on the 
comb structure, where the number of transporting par- 
ticles is not conserved, and the continuous time random 
walk (CTRW) with conservation of the total number of 
particles (transporting or not). This case of s > 1 differs 
essentially from those with s < 1 where the transport- 
ing particle move in the direction of the convection cur- 
rent j x X]. Unlike our previous consideration in the 
present study we observed analytically that an asymp- 
totic solution for the transporting particles corresponds 
to superdiffusion in the direction opposite to the convec- 
tion flow current —j x . This is a result of an irreversible 
relaxation process .9] which is "inevitable" process in the 
diffusion transport phenomena. 



First, we will describe the superdiffusion on the comb 
structure due to the inhomogeneous convection described 
by the 2D distribution function G = G(t,x,y) and the 
current 



j= { cir.<nG-LHy)^-,-D^ 
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where D5(y) and D are the diffusion coefficients for the 
x and y directions correspondingly, while the inhomoge- 
neous convection velocity is v{x,y) — v\x\ s 8{y). 

The comb model is known as a toy model for porous 
media used for the exploration of low dimensional per- 
colation clusters 01 an d an electrophoresis process 0] . 
For v = subdiffusion have been observed [IJ. A spe- 
cial transport behavior on the comb structure is that the 
displacement in the x-direction is possible only along the 
structure £-axis, say at y = 0, according to Eq. Q). Both 
the diffusion coefficient and the convection flow arc highly 
inhomogeneous in the y-dircction. There is also diffusion 
in the y-direction with a constant diffusion coefficient D. 
Therefore the Liouville equation 

- + d lvJ = 

corresponds to the following Fokker-Planck equation 



dG 
dt 



— + L FP (x)GS(y) - D 



d 2 G 
dy 2 



(2) 



with the Fokker-Planck operator of the form 



L FP (x)G 



dG 



— + sv\x\ s 1 sgn(x)G. 
ox 



The initial condition is G(0,x,y) — S(x)6(y), and 
the boundary conditions on the infinities have the 
form G(£,±oo,±oo) = and the same for the first 
derivatives with respect to x and y G' x (t, ±oo, ±oo) = 
G'yit, ±00, ±oo) = 0. The function sgn(cc) equals to 1 
for x > and —1 for the opposite case. The transport 
of particles along the structure x-axis is described by 
the function G(t,x,y = 0) = g(x,t). It should be un- 
derlined that the tails of the distributions are the most 
interesting for applications. Therefore, we are study- 
ing here the large scale asymptotics when |x| 3> 1. In 
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the Laplace-Fourier (CT) space (p, k) Eq. J5J is trans- 
formed to the following fractional equation for the func- 
tion g = g(p, k) = £F[g(t, x)] 



Dk 2 g(p, k) + ik 



d s g(p, -k) 
d\k\ s 



2^Dpg(p,k) = l. (3) 



Here the fractional Reisz derivative is the result of the 
Fourier transform |9Ul3| 



d s g(--.,-fc) 

d\k\ s 



r[\x\'9(..;X)]- 



The large scale asymptotics \x\ 3> 1 corresponds to \ k\ <C 
1 in the Fourier space. Therefore, the first term in (J3J 
can be omitted at the condition 



k 2 

hm4- = 0. 

fc->o D I 



(4) 



This approximation depends on the form of singularity 
of the the convection velocity in the limit x — > oo. It 
means that the asymptotic solutions of the homogeneous 
part of Eq © for \x\ ^> 1 depends on the exponent s in 
the power law |a;| s [l^lT^. After performing the inverse 
Fourier transform, one obtains the asymptotic of x S> 1 
solution that corresponds to the homogeneous part of Eq. 
©. It reads 



g(p> x ) = T~j7 exp 



s sgn (x) 



l)v 



(5) 



This solution describes asymptotic transport of any ini- 
tial profile. To obtain the time-dependent solution 
one carries out the inverse Laplace transform g(t, x) — 
C~ l [g(p, x)). The necessary condition theorem needs the 
negative sign of the function in the exponential in Eq. 
(jSJ. It depends only on s and the sign of v. When s < 1, 
the initial profile of particles moves in the directions of 
the convection flow, namely, in the direction of v = \v\ 
for x > 0, and v = —\v\ for x < 0. It is usual superdif- 
fusive acceleration of diffusion due to the inhomogeneous 
convection. This case together with s = 1 was considered 
in detail in 

When s > 1 the situation is much interesting and leads 
to the absolutely new effect. Indeed, for s > 1, the neces- 
sary condition to perform the inverse Laplace transform 
is v = — \v\, when x > and v — \v\, when x < 0. Hence, 
the inverse Laplace transform gives 



g(t,x) 



-sg 



n(x)D 1 / 2 \x\ 



l-2.s 



v(s — 1)V 7rf 3 



• exp 



Dx 



2-2s 



v 2 {s-l)H 



(6) 

When |x| ^> 1 and t is large enough to put the expo- 
nential to the unite, one obtains the distribution for su- 
perdiffusion of particles 



g(t,x) oc 



1 



(7) 



All moments of x higher than 2s — 2 are equal to the 
infinity. It also should be underlined that the flux on 
the infinities is vanishing. The important feature of this 
superdiffusion is that it occurs in the direction opposite 
to the inhomogeneous convection current. This new phe- 
nomenon is related to the relaxation process or it is due 
to diffusion, where the Kolmogorov conditions (see 0) 
are necessary for the inferring of the Fokker-Plank equa- 
tion (FPE). In the absence of the convection the solution 
of the FPE gives that at any moment t > the particles 
are spread over the all x-axis from the minus infinity to 
the plus infinity with the exponentially small tails. It is 
correct not only for the normal diffusion but for the sub- 
diffusive relaxation on the comb structure, as well uM 



g(t, x) 



exp 



ADu 



Du 2 



This behavior is dominate for small x even in the presence 
of the inhomogeneous convection. But for the asymptot- 
ically large x the inhomogeneous convection in the di- 
rection opposite to the spreading of particles changes the 
shape of the tail of the packet from the exponential to the 
power law in according with Eq. J7J). We call this solu- 
tion the negative superdiffusion solution or the negative 
superdiffusion approximation (NSA). 

The total number of transporting particles on the 
structure axis decreases with time 



(G) = g(t,x)dx = (4s-3)/\UDt. 



(8) 



Therefore, the distribution function JBJ describes the 
NSA when the number of particles (G) is not conserved. 
The formulation of the NSA problem with conservation of 
the total number of particles is equivalent to the case with 
a continuous distribution of the delay times where 
the total number of particles is described by the function 
Gi(t,x) = J_ oc G(t,x,y)dy. It is simply to show from 
Eq. © that 



G(t,x,y) = C 1 g(p,x)e 



(9) 



Taking this into account, one obtains the equation for G\ 
by integrating Eq. with respect to the y variable. It 
reads in the Laplace space for G\(p, x) = £{Gi(t, x)]: 

pGi(j>, x) + L F p{x)g{p, x) — 8(x) (10) 
It is simply to see from Eq. or Eq. 10 that 

L FP g(p, x) =S(x)- 2^fpDg(p 1 x) . 
Substituting this in Eq. (|10|) . one obtains that 

()(p>x) = ]^y/pDGi{p,x) . 

After substitution of this relation in Eq. (fTUfl , the CTRW 
equation in the Laplace space is 

p 1/2 G 1 + ^L FP ^)Gi=p- 1/2 S(x). (11) 
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We introduce the Riemann-Liouville fractional deriva- 
tives (see, for example, 0,0) by means of the 
Laplace transform (0 < a < 1): 



(12) 



that also implies d a [l]/dt a = 16]. Using this defini- 
tion, we write down the CTRW equation which corre- 
sponds to the comb model desribed by Eq. 



1 

2D 



(x)Gx = 0. 



(13) 



Here the initial condition is G±(0,x) — S(x). For the 
asymptotically large scale x 1 ( or i < —1), we ne- 
glect the inhomogeneous term together with the second 
derivatives with respect to x in Eq. (|llf> to obtain the 
following equation 

v(d\x\ s Gi/dx) + 2iV /2 Gi = 

with the NSA related to the CTRW by 



exp 



2Dp 1/2 £ 1 ~ s sgn (x) 



v{s - 1) 



(14) 



In the rest of the Letter we infer the NSA in the frame- 
work of the Liouville-Green asymptotic solution for lin- 
ear differential equations 0|. We show that the per- 
formed approximation due to the condition is satisfac- 
tory good and corresponds to the Liouville-Green (LG) 
approximation also called the WKB approximation [l7| . 
The CTRW equation (|13|) in the generalized form reads 



(15) 



where < a < 1. Hence, for x 1, we obtain the 
homogeneous part (lhs) of Eq. I|ll|) , where the item p 1 / 2 
is substituted by p a . It reads 

- G'l + Zx s G\ + ^V^Gi + p a G 1 = . (16) 
D D 

The term in the first derivative is removed from the equa- 
tion by the substitution 



d = exp[vx s+1 /2D{s + 1)] 



(17) 



Thus w" = R(x)w, where 



R(x) 



v 2 x 2s 
AD 2 



2sD 



-s-l 



8DD 



p x l 2 x- 2s \ 



The LG approximation for w, that satisfies to the ac- 
cepted boundary conditions (see Eq. @), is 



w 



BR- 1 / 4 exp [-jR 1 / 2 dx] 



2D(s+l) v(s-l) J ' ( 18 ) 

where B is a constant. Analogously, we obtain the LG 
solution for the negative x <C —1. Therefore, taking 

B = and a = 1/2, wc obtain that Eq. (O co- 

incides exactly with the solution 114fl . It means that 
removing the second derivatives from Lpp or, the same, 
the term k 2 in the Fourier space in the limit k — *■ cor- 
responds to the Liouville-Green approximation for the 
Fokker-Planck equation with inhomogeneous (superdif- 
fusive) convection. This asymptotic solution is superdif- 
fusive transport of particles in the direction opposite to 
the convection current, namely it is the NSA. 
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